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Fragmentation of Nuclei at Intermediate and 
High Energies in Modified Cascade Model 

G. Musulmanbekov, A. Al-Haidary 

Abstract 



The process of nuclear multifragmentation has been implemented, 
I/"") ■ together with evaporation and fission channels of the disintegration of 

excited remnants in nucleus-nucleus collisions using percolation theory 
ly-j \ and the intranuclear cascade model. Colliding nuclei are treated as 

face-centered-cubic lattices with nucleons occupying the nodes of the 
lattice. The site-bond percolation model is used. The code can be 
| applied for calculation of the fragmentation of nuclei in spallation and 

' multifragmentation reactions. 



"o ■ 1 Introduction 

The intranuclear cascade model is one of the basic tools for analyzing spalla- 
tion and multifragmentation processes in nuclear collisions. In the traditional 
^ ■ cascade model of hadron-nucleus and nucleus-nucleus interactions, particle 

production is treated in two stages. In the first fast stage, intranuclear cas- 
cade occurs inside the target and (or) the projectile nuclei and some nucleons 
from the target and projectile nuclei are knocked out, together with mesons. 
In the second stage, residual nuclei (generally, in an excited state) divide into 
two remnants in the fission channel or evaporate protons, neutrons and/or 
light nuclei, including helium isotopes. However, experimental data indicate 
that, at intermediate energies, a third competing process, multifragmenta- 
tion, comes into play, in which excited remnants break up into intermediate 
mass fragments (IMF). There are two approaches for theoretical description 
of multifragmentation: dynamical and statistical. In statistical multifrag- 
mentation models, an excited remnant achieves thermal equilibrium state 
and then expands, eventually reaching the freeze-out volume. At this point 
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it fragments into neutrons, light charged particles and IMFs. In dynamical 
models IMFs are formed at the fast stage of nuclear collision via dynamical 
forces between nucleons during the evolution of the total system of interacting 
projectile and target. In this case the whole system and its parts (projectile 
and target remnants) never pass through states of thermal equilibrium. 

There is one more approach for description of the process of multifrag- 
mentation: percolation theory. Percolation models treat the nucleus as a 
lattice with nucleons located at nodes of the lattice. It has been found that 
results of percolation calculations depend significantly upon the details of the 
lattice structure. For reasons of computational convenience, the simple cubic 
lattice has been most frequently used in multifragmentation simulations 0, 
but several studies have found P, |3| that the face-centered-cubic lattice more 
accurately reproduces the experimental distributions of fragment masses and 
their energy spectra. 

Although lattice simulations have been found to reproduce multifragmen- 
tation data surprisingly well, there has been little examination of the role 
of the lattice arrangement of nucleons inside nuclei. That is, lattices were 
employed more as computational techniques, rather than as formal nuclear 
models. The appearance of solid state models of nuclear structure can be 
dated from the paper of Pauling in 1965 @, [|, ||. The most attractive lat- 
tice model is the face-centered-cubic (FCC) model proposed by Cook and 
Dallacasa0 because it brings together shell, liquid-drop and cluster char- 
acteristics, as found in the conventional models, within a single theoretical 
framework. Unique among the lattice models, the FCC reproduces the entire 
sequence of allowed nucleon states as found in the shell model. 

In the present paper we further develop the modified intranuclear cascade- 
evaporation code, (MCAS) elaborated by one of the authors 0, with the aim 
of inclusion of multifragmentation channels. The word "modified" relates to 
the implementation of the concept of "formation time" into traditional cas- 
cade calculations, as described in Section 3. The goodness of fit of the MCAS 
to experimental data has been reported in previous papers J7|, || concerning 
multiparticle production in nucleus-nucleus collisions at intermediate and 
moderately high energies (up to 10-20 GeV/n). Since traditional cascade 
models consider nuclear structure as a dilute fermi gas, we reconstructed the 
model in the framework of the lattice nuclear model. For this purpose we 
implemented the FCC lattice arrangement of nucleons for the colliding nuclei 
according to the algorithm proposed in reference (Section 2). Calculations 
of multifragmentation channels are performed on the basis of the bond-site 
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model of percolation theory (Section 4). Comparisons with experimental 
data are given in Section 5. 

2 FCC lattice model of nuclear structure 

The FCC packing of nucleons, with protons and neutrons occupying lattice 
sites in alternating layers, can be seen as consisting of four interpenetrating 
cubes. A nearest-neighbor distance of about 2.0262 fm reproduces the known 
core density of nuclei (0.17 nucleons/fm 3 ). The essence of the geometry of 
the FCC model can be shown using the quantum numbers that are assigned 
to each nucleon in the conventional shell model M. It is known that a 
nucleon's distance from the center of the nucleus determines its principal 
quantum number n. The distance of the nucleon from the "nuclear spin 
axis" determines its total angular momentum (quantum number j). Finally, 
the distance of each nucleon from the y — z plane determines its magnetic 
quantum number m. The inherent simplicity of the FCC model is evident 
in the FCC definitions of the eigenvalues: 



n = (\x\ + \y\ + \z\ - 3)/2, (1) 
j = (\x\ + \y\-l)/2, (2) 
|m| = \x\/2, (3) 

where the sign of the m value is determined by the intrinsic spin orientation 
of the nucleon in the antiferromagnetic lattice (spin up =| and spin down 
=— ~). Conversely, the coordinate values can be determined solely from the 
nucleon eigenvalues: 

x = |2m|(-l) m+ ^, (4) 
y = ( 2 j + l-\x\)(-l) i+ i +m+1 2, (5) 
z = (2n + 3- |ac| - |y|)(-l) <+n " J " _1 , (6) 

where i is the isospin quantum number. Therefore, knowing the full set of 
eigenvalues for a given set of nucleons, the configuration of those nucleons in 
3-D space relative to the nuclear center can be determined unambiguously. 
Using the fermi coordinates of each nucleon, the mean radius of the nucleus 
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with A nucleons is defined as 



R [A] — Rnudeon + ~J ^ ] T ji (7) 

where r is the Euclidean distance of each nucleon, x 2 + y 2 + z 2 , from the 

origin and R nuc ieon is the nucleon radius. The calculated charged radii for 
various nuclei are in good agreement with experiment. 



3 Intranuclear Cascade with the Nuclear Lat- 
tice Model 

Nucleon coordinates (4 - 6) for the target (projectile) nucleus are generated 
in accordance with the algorithm given in reference [jKfl . For each nuclear 
collision, lattices of target and projectile nuclei are oriented randomly in 
relation to the collision axes. This random orientation of the nuclear lattice 
in 3-D space mimics the Woods-Saxon distribution of nuclear density for 
medium and heavy nuclei. Nucleon momenta inside the nucleus, p, are 
generated uniformly in the space < |p| < pp. The bound Fermi momentum 
Pp relates to the local nucleon density as 

p F = (3n 2 ) 1/3 hp 1/3 (r). (8) 

An inelastic collision of two nuclei is an incoherent superposition of baryon- 
baryon, meson-baryon and meson-meson elastic and inelastic interactions. 
Elastic and inelastic cross sections and kinematical features of the elastic 
scattering are taken from experiments. Description of the inelastic event gen- 
erator is given in Appendix. All interactions are arranged into four groups. 
Group C - interactions of the nucleons of the projectile nucleus with those 
from the target nucleus. All secondary particles produced in any group of 
interactions are considered as cascade particles. 

Group A - interactions of the cascade particles with the nucleons of the tar- 
get nucleus; 

Group B - interactions of the cascade particles with the nucleons of the in- 
cident nucleus; 

Group D - so called "cascade-cascade" interactions — interactions of cascade 
particles with each other. 



4 



The probability of interaction of particles % and j is defined by a black 
disk approximation: 

p(bl) = e(bl-a tot M, (9) 

where bij is the impact parameter between hadrons % and j, and o to t is their 
total cross section. Cross sections of resonances in subsequent interactions 
are taken to be the same as for stable particles. The evolution of the inter- 
acting system is considered as follows. At some instant of time t all possible 
interacting pairs in each group (A,B,C,D) are determined. Among all possi- 
ble interactions that one is chosen to be the first if it occurs before others, i.e., 
At = min{ti}; then the positions of both nuclei and all cascade particles are 
moved to new positions corresponding to a new instant of time U — > U + At. 
Since the formation of hadronic states of secondary particles takes some time, 
we apply the concept of formation time (zone) for consideration of their sub- 
sequent interactions. The formation time relates to the development of the 
cross section of the produced particle during its propagation inside the nu- 
clear medium. We use the exponential form of the evolution of cross sections 
until the subsequent collision occurs 

a l 2 = a[{\ - (1 - xV r/(7To) ) (10) 

for the leading particle, 

o™ = d?- « - x m a{)e- Tl/ ^ To) (11) 

for the m-th produced particle, where o\ is the normal cross section of the 
incident particle in the first collision, a\ is the cross section of the remnant of 
the projectile (leading particle) in the second collision, a™ is the cross section 
of the m-th produced particle, a™ is the cross section for this type of particle 
in the normal state, 7 is the Lorentz- factor and To is an adjustable parameter 
corresponding to the mean value of the formation time in the rest frame of 
the particle. For r + 1-th inelastic rescattering of the incident particle, the 
cross section is defined as 

^ + i = ^fl(l-(l-^)e- ri/(7iro) ) (12) 

i=i 

Among secondaries, s-wave resonances (deltas, rho and omega mesons) 
can be produced. The hadronic event generator is briefly described in the 
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Appendix. During the evolution of the system, the produced resonances 
may decay before their subsequent interactions. A check is made whether 
the Pauli principle is satisfied both for all interactions and for the decay of 
resonances. The cascade stage of particle generation is completed when all 
cascade particles have left both nuclei or have been partly absorbed by them. 
In this way, the first fast stage of multiparticle production of the nuclear 
collision has been completed. After replacing the fermi-gas nuclear model 
by the FCC lattice, we compared the results of simulations using both models 
on multiparticle production in intermediate and high energy nuclear collisions 
and have found that they are identical. The first measurable characteristic 
of nucleus-nucleus collisions is the reaction cross section. In intranuclear 
cascade models the cross section is defined by the ratio of the number of 
realized inelastic collisions, N in , to the total number of trials, N tnal : 

ivrin 

° reac = j^^A + R B + Ay, (is) 

where Ra,b is the radii of the colliding nuclei and A is the radius of the 
strong interaction. The fermi-gas and FCC lattice models are in agreement 
to within an accuracy of 5 percent. 



4 Fragmentation of Excited Remnants 

The number and total charge of the remaining nucleons in each remnant spec- 
ify the mass and charge numbers of the residual nuclei. In general, remnants 
are in excited states and possess angular momentum. The excitation energy 
of each remnant nucleus is determined by the energy of the absorbed parti- 
cles and the "holes" remaining after nucleons have been knocked out during 
the intranuclear cascade process. The momentum and angular momentum 
of the residual nucleus are evaluated in light of the conservation of momen- 
tum and sequentially followed for each intranuclear interaction. Thus, there 
are three competing processes for the disintegration of the excited remnant 
nucleus: evaporation, fission and multifragmentation. In the standard in- 
tranuclear cascade model, only the first two processes are taken into account 
TT| . The purpose of the present study was to implement multifragmentation 



on the basis of percolation theory and to determine the relative weights of 
the above three competing processes. We have done this by applying the site- 
bond percolation model. We assume that nucleons occupying lattice sites are 
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connected with their neighbors via bonds which schematically represent two 
body nuclear forces. In the fast stage, during the development of intranuclear 
cascade some nucleons occupying the sites of the FCC lattice of the target 
(projectile) nucleus are knocked out, leaving "holes" at those sites. We say 
that these sites are broken. The ratio of the number of broken sites to the 
total number of sites (the mass number of the target or projectile) charac- 
terizes the degree of destruction of the target (projectile) nucleus after the 
cascade stage. This ratio depends on the collision energy, the mass numbers 
of colliding nuclei and, particularly, on the impact parameter of the collision. 
In peripheral collisions, mainly peripheral nucleons are knocked out, mean- 
ing that, with high probability, the remaining nucleons form one cluster in 
which all sites are occupied. In collisions with more centrality, corresponding 
to intermediate or small impact parameters, nucleons are knocked out mainly 
from the nuclear interior and the target (projectile) remnant represents the 
lattice with some sites broken. As mentioned above, remnants, in general, are 
in excited states. The larger the impact parameter, the smaller is the num- 
ber of broken sites and the less is the excitation energy of the remnant. This 
initial condition is preferable for equilibration and thermalization of excited 
nuclear media and allows one to use evaporation and fission mechanisms for 
subsequent disintegration of the excited remnant. With increasing central- 
ity of the collision, the number of broken sites increases (large destruction), 
leading to increasing excitation energy of the remnant nucleus. For this case, 
there is no conventional understanding of the mechanism of disintegration of 
an excited remnant (thermal break-up with statistical multifragmentation, 
liquid-gas phase transition, sequential evaporation, cold shattering break-up, 
etc.). However, it is obvious that when there is considerable destruction of 
the remnant, there is no possibility for equilibration and thermalization over 
the whole volume of the remnant. 

For excited remnant disintegration we specify the bond breaking probabil- 
ity, as an input parameter pbond, in the form of impact parameter dependence: 



where Ra and R B are the radii of the colliding nuclei. This anzatz can be 
derived from considerations of the collision geometry. The cluster counting 
algorithm, developed by authors, looks for clusters (fragments): whether 
neighboring nucleons are connected via bonds or not. Only first -nearest and 




(14) 
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second-nearest neighbors are taken into account in the counting algorithm. 
In the initial FCC lattice each nucleon has 12 first-nearest neighbors at a 
distance of 2.0262 fm and 6 second-nearest neighbors at 2.8655 fm. As a 
result of this counting algorithm, we obtain the mass and charge distribution 
of the fragments. Although this approach is statistical and the probability 
of any bond to be broken does not depend on its position, the probability of 
disintegration of the remnant on multiple clusters (fragments) will be higher 
in the vicinity of regions with many broken sites. From this it follows that the 
process of multifragmentation is influenced by the dynamics of the collision, 
i.e., according to our scenario it is not purely a statistical process. 

Next we specify the energetic characteristics of the radiated fragments. 
In general, in its proper frame the remnant possesses rotational energy,i? ro *, 
and excitation energy, E*, which are used in the summation of the rotational, 
Eff, and kinetic energies of fthe ragments, E kl r n , their excitation energies, 
Ef r , and the energy of the coulomb interactions of the fragments, E c ™ 1 : 

E r °t _j_ E * E r °^ -\- J^ ia -)- Ji]* -)- E cou ^ 

£ E rot + E kin (Ai, Zi ) + E*(A t , Zi ) + W ( 15 ) 

z r t j 

In the standard intranuclear cascade model the contribution of rotational 
energy, E^f, is small compared with other terms, at least for light nuclei, 
as projectiles. Whether or not this is the case in reality is unknown. Large 
rotational energies could be realized in this approach if we included nuclear 
viscosity. In the current calculations we neglect the first term. Moreover, for 
computational convenience we make additional simplifications in Eq. (15). 
Since the coulomb repulsion of the charged fragments increases their kinetic 
energies, we define the resulting kinetic energies of the fragments as follows: 

E Hn + E coul = E kin (A h Zl ) + l-J2— = Yl £ ^ Z *)' ( 16 ) 

Z Tij 

Another simplification concerns the excitation of the fragments: we assume 
that only one fragment among others is excited, the mass number of which 
is maximal. It is justified, particularly, when comparing model with data 
obtained through inverse kinematics because the experimental setup registers 
a majority of radioactive fragments as well. Therefore, the excitation energy 
of the remnant, E*, is converted into the kinetic energies of the fragments 
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and the excitation energy of the fragment with maximal mass: 



E* 



E*(A 



max j ^max) 



(17) 



With these simplifications we generate the energy distribution of frag- 
Before the collision 



ments applying considerations proposed in reference [12 



the nucleons have a momentum distribution that is uniform inside the Fermi 
sphere of radius pp. After the collision the distribution in the vicinity of the 
beam propagation is wider because of intranuclear interactions accompanied 
by local excitation of the nuclear medium. This can be written in the form: 



n(e) oc 1/ < 1 + exp 



e — e F 



T, 



eff 



(18) 



where e = p 2 /2m, and e F is the boundary Fermi energy. The "effective 
temperature" is given by 



T. 



eff 



cE*/N x 



br ■ 



(19) 



where c is an adjustable parameter and Nf, r is the number of broken sites. 
Kinetic energies of nucleons composing the fragment are generated according 
to distribution (18), and, summing up all vector momenta directed randomly 
in 3-D space, we obtain the momentum of the fragment. In such a way, we 
generate momenta of all produced fragments. The remaining part of the 
remnant excitation energy, (17), is assigned to the fragment with maximal 
mass number. And, of course, we take into account the conservation of energy 
and momentum for the whole reaction. 



5 Comparison with Experiment 

Observation of residues emerging from spallation reactions in direct kine- 
matics still remains a difficult task. Collisions of protons and light nuclei 
with heavy ions performed at GSI in inverse kinematics allows one to de- 
termine the production of residues prior to j3 decay. This provides a good 
opportunity to compare the available data with theoretical models to achieve 
a better understanding of the mechanisms of reactions which is, today, far 
from satisfactory. Until now calculations have been performed by different 
versions of intranuclear cascade followed by the evaporation model. As seen 
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Figure 1: Mass distribution of residues produced in 0.5 GeV proton induced 
197 Au; Pbond = 0.5. Data are from paper [[ITJ . 
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Figure 2: Mass distribution of residues in the reaction 197 Au + p at 800 
A MeV; pbond = 0.57. Dashed curve is outcome of calculations without the 



contribution of multifragmentation channels . Data are taken from paper 1 14 . 
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Figure 3: Mass distribution of residues in the reaction 208 Pb + p at 1 A GeV; 
Pbond = 0.6. Data are taken from paperjl5 . 
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Figure 4: Mass distribution of residues in the reaction 208 Pb + d at 1 A GeV; 
Pbond — 0.62. Data are taken from paper[|TB]. 
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Figure 6: Isotopic distribution of spallation residues in reaction 208 Pb + p at 
1 A GeV; Pbond = 0.6. Data are from paper [|T5|. 
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from the previous section, our model also includes multifragmentation chan- 
nels in the framework of the percolation approach. Here we define the values 
of the input parameters, which are the bond breaking probability, Pbond{fy, 
in Eq. (14) and the constant c in Eq. (19). Pbond{fy depends on the energy 
of the collision and the type of reaction. For a specific reaction, it is obvi- 
ous that site and bond breaking probabilities are small at low energies and 
start growing with increasing energy, reaching constant values at the regime 
called "limiting fragmentation". Limiting fragmentation is reached at differ- 
ent energies for different reactions. Therefore, at low energies the dominating 
mechanisms of disintegration of excited remnants are evaporation and fission. 
As the energy of the collision grows, the contribution of multifragmentation 
processes increases, depending on the site and bond breaking probabilities. 
Since the number of broken sites is defined automatically during the develop- 
ment of the intranuclear cascade, only the bond breaking probability remains 
to be input as a parameter. For proton induced reactions, p& on d(0) changes 
from 0, at an incident proton energy of a few tens of MeV, to 0.77 at the 
limiting fragmentation energy (3-4 GeV). With regard to the parameter c 
in Eq. (19), its value is chosen to be 0.7, energy independent, for all types 
of reactions. Comparison of the model calculations for mass distributions in 
spallation reactions pAu at 0.5 and 0.8 GeV, pPb at 1 GeV, and dPb at 2 
GeV are shown in Figs. 1-4. As seen from the figures, at energies lower 
than those corresponding to the limiting fragmentation regime, mass yield 
distributions of residues in proton-induced reactions have well-pronounced 
bell-shaped curves in the central part, corresponding to the contribution of 
fission channels. Evaporation channels give a dominating contribution in the 
right peak of the distribution with a plateau at high mass residues. The 
values of the level density parameters for evaporation and fission are taken 
to be 0.1A MeV" 1 , the same for both and independent of the type of re- 
action and collision energy. Figures 5-7 show the isotopic distributions of 
residues produced in the following reactions: pAu at 0.8 GeV, pPb at 1 GeV 
and deuteron - Pb at 2 GeV. Calculated distributions are shifted toward 
neutron-rich isotopes for lighter residues. We think that underestimation 
of the proton-rich isotopes and overestimation of the neutron-rich isotopes 
could be corrected by better description of proton-neutron competition in 
fission-evaporation channels. 

As the bombarding energy increases, so does the contribution of multi- 
fragmentation channels and, correspondingly, the contribution of evaporation 
and fission channels decreases. This leads to filling of the dips at both sides 
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Figure 8: Mass distribution of fragments in the reaction p + 108 Ag at 11.5 
GeV; pbond = 0.77. Data are taken from paper |T7 . 
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Figure 9: Mass distribution of fragments in the reaction p + 197 Au at 11.5 
GeV; pbond = 0.77. Data are taken from paper [1181. 
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Figure 10: Mass distribution of fragments in the reaction 12 C + 108 Ag at 4.5 
A GeV; Pbond — 0.72. Data are taken from paper [[01 . 



of central bell-shaped curve and a decrease in the height of the right peak. 
As mentioned above, rising collision energy leads to an increasing number 
of broken sites and bonds in the nuclear lattices. This, in turn, results in 
increasing yield in multifragmentation channels. This tendency is already 
evident in reaction Au + p at 0.8 A GeV (Fig. 2). The mass distributions of 
fragments in proton- and carbon-induced reactions on Ag and Au at energies 
corresponding to the limiting fragmentation regime are shown in Figures 8 - 
10. 



6 Conclusions 

A new version of the Modified Cascade Model for intermediate and high 
energy nucleus-nucleus collisions including multifragmentation channels has 
been developed. Colliding nuclei are represented as face-centered- cubic 
lattices. Multifragmentation is calculated in the framework of percolation 
theory with usage of a site-bond percolation model. This version is able to 
reproduce reasonably well both spallation and multifragmentation processes. 
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7 Appendix: Hadronic Event Generator 



Monte Carlo simulation of inelastic events is performed in several steps. The 
first step in the generation of an exclusive event is the evaluation of the initial 
cm. energy portion available for production of secondaries 

W = ^Ei = k^Ts (20) 

where Ei is the energy of the i-th particle (excluding leading particles), k 
is inelasticity. Fluctuation of the inelasticity from event to event leads to 
the distribution P(k). There are no comparable theoretical methods for 
calculation of P{k). It has been shown in Ref. |20] that one may fit the 
inelasticity distribution with a beta distribution 

P{k, s) = k a -\l - kf- 1 /B{a, b) (21) 



B(a,b) = T(a)T(b)/T(a,b) (22) 

< k(s) >= a/(a + b) (23) 

where T(a), T(b) and T(a,b) are gamma functions; s is the dependence of 
P(k, s) and < k(s) > is enclosed in parameters a and b. Up to the ISR 
energies one can neglect this s-dependence. In the second step the energy 
W is distributed between secondary particles whose kinematical character- 
istics are generated in correspondence to a cylindrical phase space model. 
Parameters of the cylindrical phase space model are adjusted by comparing 
the results of simulation of pion-nucleon and nucleon-nucleon interactions 
with experimental data. The remaining part of cm. energy (1 — k)y/s is 
distributed between remnants of the interacting particles (so-called leading 
particles) according to the conservation of energy and momentum. 

Pi+Pn = J2 Ti ( 24 ) 



Ei + E n = (1 - k)y/a (25) 

where Pi is the momentum of the i-th produced particle, and Pj, Pn and 
Ei, En, are momenta and energies of leading particles. Interacting nucleons 
(mesons) can transform into nucleons (mesons) and s - wave resonances (A 
- isobars and p, u - mesons). Transition probabilities are calculated with 
the use of a one-pion-exchange (OPE) model. 
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